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; ABSTRACT 

I derive a nonlinear local relation between the redshift-space density field and the real- 
, ^ \ space velocity field. The relation accounts for radial character of redshift distortions, 

and it is not restricted to the limit of the distant observer. Direct comparisons between 
the observed redshift-space density fields and the real-space velocity fields possess all 
of the advantages of the conventional redshift-space analyses, while at the same time 
they are free of their disadvantages. In particular, neither the model-dependent re- 

■ construction of the density field in real space is necessary, nor is the reconstruction 
^ \ of the nonlinear velocity field in redshift space, questionable because of its vorticity 

at second order. The nonlinear redshift-space velocity field is irrotational only in the 
distant observer limit, and that limit is not a good approximation for shallow catalogs 
QQ ■ of peculiar velocities currently available. Unlike the conventional redshift-space com- 

' parisons, the comparison proposed here does not have to be restricted to the linear 

regime. Accounting for nonlinear effects removes one of the sources of bias in the esti- 
I mate of (3. Moreover, the nonlinear effects break the Omega-bias degeneracy plaguing 

■ all analyses based on linear theory. 
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C3 1 INTRODUCTION 
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_ Comparisons between density fields of galaxies and the fields of their peculiar velocities are a powerful tool to measure the 

cosmological parameter Q.. This is because in the gravitational instability paradigm the density and the velocity fields are 
$_( tightly related, and the relation between them depends on fi. 
. 5^ , These comparisons are commonly performed in real space. A necessary ingredient of real-space analyses is the recon- 

struction of the galaxy density field in real space from the observed galaxy field, which is given only in redshift space. While 
performing such a reconstruction, one has to correct the redshifts for peculiar velocities of galaxies. However, the amplitude 
of these velocities depends on f2, or in the case of bias, on /3. This is a serious drawback of real-space comparisons: they 
necessarily assume the value of a parameter which is to be subsequently estimated. 

To avoid this problem, the comparisons in redshift space have been proposed (Nusser & Davis 1994, hereafter ND). 
However, the velocity field in redshift space is irrotational only at the linear order,^ so that the redshift-space analyses 
must be restricted to the linear regime. This is unsatisfactory, since the derived amplitude of the density fluctuations from 
current redshift surveys (e.g.. Fisher et al. 1995) and from the potent reconstruction of density fields (Dekel et al. 1990 
and Bertschinger et al. 1990), slightly exceeds the linear regime. For example, the density contrast in regions like the Great 
Attractor or Perseus-Pisces is about unity even when smoothed over scales of 1200 km s~^ (Sigad et al. 1998). Future redshift 
surveys and peculiar velocity catalogs are expected to provide reliable estimates of density and velocity fields on scales 
where nonlinear effects are certainly non-negligible. Those eff'ects may lead to interesting consequences, such as breaking the 
degeneracy between f2 and bias (Dekel et al. 1993, Bernardeau et al. 1999, Chodorowski 1999; hereafter C99), and therefore 
they should be accounted for. 



* Chodorowski & Nusser (1999) have shown that the nonlinear redshift-space velocity field is irrotational in the distant observer limit. 
However, the catalogs of peculiar velocities are not yet deep enough for this assumption to be well satisfied (Zaroubi &; Hoffman 1996, 
Szalay, Matsubara & Landy 1998). 



© 0000 RAS 



2 M. J. Chodorowski 



Here I propose to compare the redshift- space density field directly to the real-space velocity field. Such a comparison 
enables one to avoid the model-dependent reconstruction of the density field in real space. Also, the vorticity of the velocity 
field is no longer a problem, because (before shell crossings) the real-space velocity field is irrotational without any restrictions. 
Therefore, the comparison proposed here combines all advantages of real-space and redshift-space analyses, while at the same 
time it is free from their disadvantages. 

The present analysis is based on the work of C99, who derived a similar relation in the distant observer limit (DOL). 
However, as pointed out first by Zaroubi & Hoffman (1996), and more recently by Szalay, Matsubara & Landy (1998), the 
DOL approximation is not well satisfied by available catalogs of galaxy redshifts and peculiar velocities. Here I relax this 
assumption, and I take fully into account the radial character of redshift distortions. I show how the derived relation can 
be inverted, i.e., how to reconstruct the real-space velocity field from the redshift-space density field. This is a necessary 
ingredient of the so-called velocity-velocity comparisons. The results of the paper are then compared to the results of ND. 

The outline of this paper is as follows: in Section ^ I derive the exact relation between the redshift-space density and 
the real-space velocity. In Section ^ I expand it explicitly up to linear and second-order terms. Next, in Section ^ I derive a 
similar relation between the galaxy density field and the velocity field under an assumption of a nonlinear but local bias. In 
Section |5| I show how such a relation can be used to break the f2-bias degeneracy in velocity-velocity comparisons. Summary 
and conclusions are in Section M. 



2 EXACT RELATION 

The transformation from the real space coordinate, r, to the redshift space coordinate, s, is (Kaiser 1987; hereafter K87) 

S = r + VrT , (1) 

where Vr{r) — v - r; v{r) is the real-space velocity field, and velocities are measured relative to the Local Group. The relation 
between the redshift-space mass density field, 5s, and the real-space velocity field can be obtained from the equation 



^ n! 

n=0 



(j)[s{r)] 



(2) 



Here, is the selection function, 5{r) is the real-space density field, and J(r) is the Jacobian of the transformation (^), 
J{r) = {l + vr/r)\l + v'r), (3) 
where ' = d/dr. 

Equation (^) follows from the continuity equation; its derivation is outlined in Appendix The equation is exact; however 
it is strictly valid only in the absence of triple-value zones in redshift space. It is a generalization of equation (A6) of C99, 
which was derived in the DOL. 

The relation between the redshift-space density and the real-space velocity follows from equation ^ when the density- 
velocity relation (DVR) in real space is used. The DVR in real space is well known, it has been extensively studied both 
analytically and numerically (Nusser et al. 1991, Bernardeau 1992, Gramann 1993, Mancinelli et al. 1994, Mancinelli & Yahil 
1995, Chodorowski 1997, Chodorowski & Lokas 1997, Chodorowski et al. 1998, Bernardeau et al. 1999, Canon et al. 2000, 
Kudlicki et al. 2000). While a number of accurate formulas has been proposed, for the purpose of this Section it is sufficient 
that the real-space density is a well-known nonlinear function of the first-order velocity derivatives, with an explicit dependence 
on SI, 

5{r)=r[dv,/dr,{r)M. (4) 
If the selection function is given by a power law, (f> (x , then 



If not, the expression (j){r) / cj>{s) should be expanded explicitly. In any case, it will be a function of Vr/r. Therefore, both 
(j>{r)/(l>{s), 5 and J are functions of the real-space velocity field and its derivatives, and equation (^ becomes a relation 
between the redshift-space density and the real-space velocity. 

The key point of equation ^ is that it relates the redshift-space density at a point s to an expression involving the 
real-space velocity field evaluated at f = s. (r is a real-space point, in general different from r, which is related to s by eq. ^.) 
This has been achieved by effectively expanding real-space quantities, appearing in the continuity equation as functions of r, 
around f = s (for details see C99 and Appendix This is very convenient, since now we can treat the two fields as if they 
were given in the same coordinate system. In other words, we can compare the redshift-space density at a redshift, say, 3000 
km s~^, directly to the real-space velocity at a real-space position 3000 km s~^. 
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Relation (Q) involves an infinite series in velocity. If the density-velocity comparison is performed on scales large enough 
so that they are only weakly nonlinear, then this series can be truncated. Truncation at order q means that the terms of the 
order g + 2 in velocity or higher are neglected. The truncation procedure should be done with care, because nonlinear effects 
are stronger in redshift space than in real space (e.g., Taylor & Hamilton 1996, Scoccimarro, Couchman & Frieman 1998, this 
paper). For example, the formalism presented here fails for triple- value zones (regions with shell crossings in redshift space). 
However, these regions or structures in real space may be just after turning around, quite a mildly nonlinear condition. Thus 
far, all density-velocity comparisons in redshift space have been restricted to linear regime (e.g., ND, Fisher et al. 1995). In 
the next section I will derive an explicit expression for the relation up to second-order terms. Moreover, it will become 
clear how to extend the derivation to still higher order, if necessary. 



3 APPROXIMATE RELATION 
3.1 Linear relation 

To derive the linear relation, we linearize the ratio of the selection functions, the real-space DVR, and the Jacobian. We have 



(r) ^ dlnfj} Vr 



(fi[s{r)] dlnr r ' 

S^-f-'Vr-V, (7) 

and 

J ~l + v'r + 2vr/r . (8) 
Hence, 

(. + i)-j.-,rV...-.;-f2 + ^)!;i, (9) 



(/)[s(r)] ■ \ dlnr 

and we see that in the series only the term with n = contributes to linear order. Therefore 

\ amr / r 



(10) 



l~DiVr/r+{Dl^D2){vr/r)\ (12) 



This equation is equivalent to equation (3.3) of K87. It also coincides with equation (6) of ND. However, the equation 
of ND describes the relation between the redshift-space density and the redshift- space velocity, u(s), and all derivatives are 
taken there with respect to s. The two equations coincide, because 

du,/ds, = dv^/dr, + O {v'^) . (11) 
Clearly, the equations would be different at the second order. 

3.2 Second-order relation 

In the series (P), only the terms with n — and n = 1 generate the terms up to quadratic in velocity. To second order, 

0[s(r)] 
where 

i^i = ^ , (13) 
dmr 

and 

X" 2 

'»> 

Next, 

J ^l + v', + 2vr/r + 2v',vr/r + {vr/rf . (15) 

Finally, the real-space density is up to second order a local function of the two velocity scalars: the expansion (the divergence), 
9, and the shear, E. Specifically (Chodorowski 1997; see also Gramann 1993, Catelan et al. 1995, Mancinelli & Yahil 1995) 

5{r) = -r'e{r) + [e\r) - |E^(r)] . (16) 

Here, 
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Yj = YjijYjij , (l*^) 
E.^ = i {dv,/drj + dvj/dr,) - 6) , (18) 

e = Vr-V, (19) 

and I have neglected the weak ^-dependence. The symbol denotes the Kronecker delta. 

From equations (^^, ([l5|), and ( [l6| ) the contribution from the term n = in the series can be obtained. As for the 
term n = 1, there is an extra factor Vr in front of the expression in round parentheses {4>{r) /(j}{s) . . .). Therefore, it is sufficient 
to adopt for this expression its linear approximation After a straightforward but rather lengthy algebra we obtain 



Ss{s) = {-r'e-v:.-i2 + D,)vr/r+[vrir'e + v',)]' + ±r^{e^~ 

+ (2 + Di) {f-'e + 2v:.) vr/r + (l + Di + - Da) {vr/rf} 



2^ 



(20) 



In the DOL, Vr/r — » and lines-of-sight are taken at a fixed direction, so Vr — > Vz and d/dr — > d/dz, and the above 
formula reduces to equation (13) of C99. 

K87 derived a relation between the redshift-space density and the real-space velocity only at the linear order. To derive a 
related relation, between the redshift-space density and the redshift-spa.ce velocity, ND applied the Zel'dovich approximation. 
However, they subsequently restricted their analysis again to the linear regime. The reason they did so was that the redshift- 
space velocity field is rotational at second order, so that it cannot be reconstructed from radial components of velocities solely. 
The analysis presented here is free of this problem since we relate the redshift-space density field directly to the real-space 
velocity field, which is irrotational without any restrictions (before shell crossings in real space). 

Applying linear theory in density-velocity comparisons leads to a bias in the inferred value of Q. We can estimate the 
amplitude of this bias using equation (^^. For simplicity, adopt the DOL and consider a spherical top-hat overdensity. Then 

i + i)r'o+(47 + i + ^)r'o' (21) 



3/ " \^21 3 9 ^ 

(cf. eq. 17 of C99). Assuming the second term on the right-hand side of the above equation to be a small correction it is 
straigthforward to show that the relative systematic error of Q, is 

n 3 V21 3 9 y ^ ' 

In the Mark III catalog, the r.m.s. value of 9, when smoothed over scales of 1200 km s-\ is ~ 0.3 (Sigad et al. 1998). Hence, 
the relative bias of the linear estimate of 57 in redshift space is about 40% for fl = 0.3, and about 30% for f2 = 1.0. This may 
be compared to the corresponding bias of the linear estimate of Q in real space, about 20 and 10% respectively. (It is obtained 
by setting the second and the third terms in parentheses in equation [ii] equal to zero.) Thus, the nonlinear effects in redshift 
space are stronger than in real space. Obviously, the bias will be respectively greater for smaller smoothing scales. 

Equation ( |20| ) can be used to reconstruct the real-space velocity field from the associated redshift-space density field. 
Since the real-space velocity field is irrotational, it can be described as a gradient of the velocity potential, 

v(r) = -Vr^^ . (23) 

Equation ( |2c| ) reduces then to a nonlinear differential equation for the velocity potential. As long as the nonlinearities in ( |2o| ) 
are weak, an iterative method of solution can be applied. First, we solve the linear part of (po|), and the method of ND (the 
decomposition in spherical harmonics) is fully applicable. Next, we find the second-order solution by solving again the linear 
equation, with the source term resulting from the density modified by nonlinear contributions approximated by first-order 
solutions. Specifically, 

/-^A.$(^) + ^ + - 5.(5 = r) - A/-2 [$«(r)] , (24) 

where 7V2 is a sum of all terms quadratic in velocity in equation (po|), expressed as functions of the potential. 

In this way, from the redshift-space density field we can reconstruct the real-space velocity potential, and hence the 
real-space velocity field. The latter can then be compared to measured radial velocities of galaxies. Note that this is a 
velocity-velocity comparison; the application of the 'redshift-real' DVR to density-density comparisons has been discussed in 
C99. 



4 GALAXY DENSITY VERSUS VELOCITY 

Galaxies are good tracers of the velocity field induced by the mass distribution. On the other hand, they are biased tracers 
of the mass distribution itself (see theoretical arguments of Kaiser 1984, Davis et al. 1985, Bardeen et al. 1986, Dekel & 
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Silk 1986, Cen & ©striker 1992, Kauffmann, Nusser & Steinmetz 1997 and observational evidence by Davis & Geller 1976, 
Dressier 1980, Giovanelli, Haynes & Chincarini 1986, Santiago & Strauss 1992, Loveday et al. 1996, Hermit et al. 1996, Guzzo 
et al. 1997, Giavalisco et al. 1998). In this section I will derive a relation between the redshift-space galaxy density field and 
the real-space velocity field under the assumption of a nonlinear but local bias. This is only a toy model for bias because there 
are good reasons to believe that bias is in fact somewhat stochastic (Dekel & Lahav 1998, Pen 1998, Tegmark & Peebles 1998, 
Blanton et al. 1998, Blanton et al. 1999, Tegmark & Bromley 1998). However, it allows for a number of important conclusions 
to be drawn. 

Equation ^ implicitly assumes no bias between the distribution of galaxies and mass, so that the densities appearing in 
it are in fact the galaxy densities Si^^ and 5'^-' (redshift- and real-space galaxy density contrasts respectively). Hence, a more 
general form of equation (H) is 



(-1)" 



n=0 



1 a" 

y 2 ^ 



2 n 

r V,. 



<t>[s{r)] 



(g) 



J 



(25) 



In a local bias model, the real-space galaxy density contrast is assumed to be in general a nonlinear function of the mass 
density contrast (Fry & Gaztafiaga 1993; see also Juszkiewicz et al. 1995), 



Given the real-space DVR 
(-1) 

n=0 



1 9" 



equation ( |25| ) yields 
Hr) 



2 n 

r V,. 



<t>[s{r)] 



[3or{dv^/drj) + 1] - J 



(26) 



(27) 



Up to second order, T is given by expression (^|) and 

B [5{r)] = b5{r) + ifoa [5\r) - a|] . (28) 

Here, b and 62 are respectively the linear and nonlinear (second-order) bias parameters. The term o"| = {5^) ensures that the 
mean value of 5'^' vanishes, as required. Combining (^8|) and (p^), and truncating the series at second-order terms we obtain^ 



[^P'^e - - (2 + D,)Vr/r + [Vrifi-'e + V'r)] ' + + ^) Z^"' 



+ (2 + Di) {r^e + 2v'r) vr/r +{l + Di + Dl- D2) {vr/rf 
Expressing velocity in terms of the velocity potential yields 

= &^^\~s = r) ~ N2 [$<^'(r),/3,/32] 



/3-iA.<I.i^' + ^ + 



Qj.2 



2 + Di 

r dr 



where 
and 

( 4 



+ 



b2_y\ 



\21b 262/ 



(29) 

(30) 
(31) 
(32) 



5 MEASURING O. AND BIAS SEPARATELY * 

Based on (|30|), we can reconstruct the real-space velocity field from the associated redshift-space galaxy density field. Com- 
parison of the predicted velocity field to the observed one will yield the best-fit values of the parameters (5 and /32 . They are a 
combination of three physical parameters: Q., h, and 62 • Therefore, we have two constraints for three parameters, so we need 
an additional constraint on them. As this constraint one can adopt the large-scale galaxy density skewness, which involves b 
and 62 (B99). 

Galaxy skewness can be measured only in redshift space. Therefore, we need a theoretical relation between the redshift- 
space galaxy density skewness S^^'' (which we can measure) and the redshift-space mass density skewness 5*33 (which we can 
compute). Approximately this relation is (see C99 for details), 

t For simplicity, I have approximated the quantity 9^ — cr^ by 8^ — (3/2)S2. This is legitimate, since the mean value of the shear scalar 
given the velocity divergence, (E^)!^, is up to second order equal to (S^) = (2/3)crg (Chodorowski 1997). 
t This Section is similar to Section 5 of C99. I include it here for the completeness of the present paper. 
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Si 




(33) 



The redshift-space mass density skewness was computed for scale- free power spectra by Hivon et al. (1995). They found 
that it depends only very weakly on Q and that it can be well approximated by the real-space skewness, except in the case 
of n approaching unity and the spectral index n J> 0. For more realistic power spectra the calculation can be done in an 
analogous way, and the result is also expected to be well approximated by the corresponding real-space value. 

The redshift-space galaxy density skewness of the IRAS density field was measured directly by Bouchet et al. (1993). 
Kim & Strauss (1998) pointed out that sparse sampling of the IRAS galaxies, coupled with the moments method used by 
Bouchet et al., makes the skewness estimates biased low. However, they proposed a method to measure the skewness by fitting 
the Edgeworth expansion to the galaxy count probability distribution function around its maximum, properly accounting for 
the shot noise. Using mock catalogs, Kim & Strauss showed that this estimate of skewness is robust to sparse sampling; they 
then used this method to measure the skewness of the IRAS density field. Hence, we know the value of the skewness of the 
distribution of the IRAS galaxies (and we know an effective means to measure it in other galaxy catalogs). 

Equations (|3ll)-(|3^) are three independent constraints for the parameters Q, b and ^2, so that they can be used to measure 
n and bias separately. The additional constraint on Q and bias (the skewness) is to be inferred from the density field alone, 
making any additional observations unnecessary. In other words, there is enough information in the density field and the 
corresponding velocity field to break the fi-bias degeneracy. 



6 SUMMARY 

I have derived a nonlinear local relation between the redshift-space galaxy density field and the real-space velocity field. A 
direct comparison of the redshift-space density directly to the real-space velocity makes the model-dependent reconstruction 
of the density field in real space unnecessary. Furthermore, the comparison need not to be restricted to the linear regime 
because the nonlinear real-space velocity field, unlike the redshift-space one, is irrotational. 

The analysis presented here is an extension of the work of C99, who performed a similar analysis in the distant observer 
limit. However, available catalogs of peculiar velocities are too shallow for this limit to be approached (Zaroubi & Hoffman 
1996, Szalay, Matsubara & Landy 1998). Here I have relaxed this approximation, and I have taken the radial nature of redshift 
distortions fully into account. C99 showed how to apply the 'redshift-real' density-velocity relation to the so-called density- 
density comparisons. Here I have shown how to apply it to velocity-velocity comparisons as well. Namely, I have derived an 
equation enabling one to reconstruct the real-space velocity field from the redshift-space galaxy density field. I have explicitly 
derived the density-velocity relation up to linear and second-order terms, but the derivation can be straightforwardly extended 
to any desired order using the exact equation (|2^). I have compared the linear relation to that of K87 and of ND and found 
them to be in agreement. 

Applying the second-order relation derived here, or its higher-order extension, to density-velocity comparisons will remove 
one of the sources of bias (e.g., the nonlinear effects) in the estimate of /3, or Q. I have shown this bias to be substantial in the 
linear theory. For the fields smoothed over scales as large as 1200 km s~^, the relative systematic error of the linear estimate 
of Q is 30-40% (depending on the value of fl). Moreover, combining the nonlinear relation with an additional measurement 
(e.g., the skewness of the galaxy density field) breaks the degeneracy between il and the galaxy bias. 

The ultimate goal of density-velocity comparisons is to estimate the value of fl separately from the galaxy bias. The 
method proposed here offers such a possibility, and it will be developped in forthcoming papers. Issues remaining to be 
addressed are the spatial scales at which a perturbative approach in redshift space can be applied, the stochastic character 
of bias, and effects of field smoothing. Kudlicki & Chodorowski (in preparation) study the first issue. Taruya & Soda (1998) 
derived a mildly nonlinear galaxy density-mass density relation in a stochastic bias model. The calculation of the galaxy 
density-velocity relation in this model will be to some extent similar. Field smoothing affects the real-space density-velocity 
relation in a very weak way (Chodorowski & Lokas 1997, Chodorowski et al. 1998, Bernardeau et al. 1999). The effects of 
smoothing are somewhat stronger in redshift space (C99). In any case, they should be studied individually for any particular 
comparison, because of different smoothing schemes used. For example, in the analysis of ND, both the density and the velocity 
fields are smoothed, while in the velmod analysis of Willick & Strauss (1998), only the IRAS density field is smoothed. In 
accounting for smoothing, the method of C99 may be applied. 
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APPENDIX A: DERIVATION OF EQUATION (§) 

Here I outline the derivation of the exact equation expressing the redshift-space density field at s in terms of the real-space 
quantities at f = s. The derivation follows to some extent the lines of Appendix A of C99, where a similar equation has been 
derived in the DOL. 

From the conservation of the number of galaxies we have 
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where if) is the selection function and J is the Jacobian of the mapping from real to redshift space. The Fourier transform of 
the redshift-space density contrast, using the mapping (hi), is 



5s{k) 



e Os{s) 



{2tt 



d'V 



ifc-r —ik^~Vr(r) 



^[s{r)] 



[S{r) + 1] - J(r) 



(27r)» 

(cf. eq. 4 of Scoccimarro et al. 1998 in the DOL). The inverse Fourier transform of the above equation is 

d^r d^fc 



%[S{r) + 1] - J(r) I e-*'-"'-('-)e*'=-(''~'-' 



(27r)3 [ (t>[s{r)] ' 

Here, I have changed the order of integration. Expanding the first exponent of the integrand we have 



<(r) d" ^,k-(s-r) 

n! dr" ' 



hence 



<P[s{r)] 



[5(r) + l]-J(r) — 



ik-{s — r) 



dr" J (27r)3 



The integral over k yields the Dirac delta distribution, 5d(s — r). Now the following lemma will be useful: 
which can be easily proven by integrating by parts. Using this lemma in equation (A5) yields equation (^) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 
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